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We extend a simple model of a charge trap coupled to a single-electron box to energy ranges and 
parameters such that it gives new insights and predictions readily observable in many experimental 
systems. We show that a single background charge is enough to give lines of differential conductance 
in the stability diagram of the quantum dot, even within undistorted Coulomb diamonds. It also 
suppresses the current near degeneracy of the impurity charge, and yields negative differential lines 
far from this degeneracy. We compare this picture to two other accepted explanations for lines in 
diamonds, based respectively on the excitation spectrum of a quantum dot and on fluctuations of 
the density-of-states in the contacts. In order to discriminate between these models we emphasize 
the specific features related to environmental charge traps. Finally we show that our model accounts 
very well for all the anomalous features observed in silicon nanowire quantum dots. 



I. INTRODUCTION 

The diversity of nanostructures available for physicists 
to perform transport experiments at low temperature has 
grown considerably in the last decade. Following pio- 
nnering works on single-electron transistors (SETs) and 
quantum dots in the early 90's, it is now very common 
to observe and analyze Coulomb blockade in novel struc- 
tures such as carbon based electronics or semiconducting 
nanowires. We have developped in the past [l| a silicon 
nanowire MOSFET which turns into a remarkably sta- 
ble and simple SET below approximately 10 K. As it is 
widely accepted that background charges in the vicinity 
of the Coulomb island or the tunnel barriers are respon- 
sible for the large 1/f noise observed in many SETs, we 
developped a model based on a single charge trap ca- 
pacitively coupled to the SET. This picture has been 
considered before in metallic single electron transistor, 
but without predictions for the Coulomb blockade spec- 
troscopy ^ . A peculiarity of our model is that the charge 
trap energy is not only sensitive to both source-drain and 
gate voltages, but also to the dot occupation number, be- 
cause of the capacitive back-action of the dot upon the 
trap. We have shown that solving the rate equations 
for this model gives the same sawtooth pattern than 
observed experimentally [3]. Sawtooth-like distortions 
of Coulomb diamonds are widely reported in the liter- 
ature of quantum dots, for instance in carbon nanotubes 
SETs [4], graphene 0, molecules in gaps Ig, M or epi- 
taxial nanowires Q. As the data published on Coulomb 
blockade transport spectroscopy, either from us or other 
groups, have released more and more details, we have 
pushed our model further in order to see if the features 
observed by experimentalists and usually attributed to 
other effects can also arise from our model. The most 
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FIG. 1: Equivalent circuit for our model of a charge trap (red) 
coupled to a quantum dot (black). The trap occupation is 
or 1 electron and it carries a very small current compared to 
the dot. The dot is treated within the orthodox model. 
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common feature is certainly lines of differential conduc- 
tance above the blocked region, parallel to the diamonds 
edges. These lines are often observed in undistorted dia- 
monds, therefore one can conclude that our Background 
Charge (BC) model cannot apply in this case, as it pre- 
dicts a clear distortion of the pattern. In this work we 
show that with the right choice of capacitive couplings 
our BC model does predict these lines in Coulomb spec- 
tra. It also gives negative differential conductance lines in 
a natural way without introducing adhoc hypothesis. We 
compare the BC model to the excitation spectrum (ES) 
model and the density of states (DOS) fluctuations model 
which are both usually invoked to explain such lines. We 
explain in details the origin of the most important fea- 
tures predicted by the BC model, the differences between 
the 3 models and how they can be distinguished. Finally 
we report data on silicon nanowire transistors where the 
traps are implanted arsenic dopants. We observe lines 
of differential conductance which are well reproduced by 
the BC model. 
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FIG. 2: Simulated Coulomb diamonds (differential drain-source conductance versus drain and gate voltages) for our model 
depicted in Figured] Cgt = 0.004 aF <^ Cgd = 16 aF, i.e. the gate-trap capacitance is very small compared to the gate-dot 
capacitance. Near trap degeneracy (Vg = V) Coulomb diamonds are replicated because of the oscillating trap occupancy with 
Vg (see Fig. [3]). At larger Vg the diamonds are shifted and undistorted; the trap is always occupied at zero bias but can still be 
empty at finite bias. This gives lines of differential conductance (see Fig. |4| slowly evolving with Vg. 



II. THE BACKGROUND CHARGE MODEL 

In this paper we consider the general case represented 
in figure [1] of a quantum dot connected to a source, a 
drain and a control gate allowing to shift its electrostatic 
energy. In addition a charge trap occupied by or 1 elec- 
tron is located nearby. The current through the trap is 
negligible compared to the main current through the dot 
but the trap is weakly coupled by tunneling barrier to one 
contact and to the dot to allow its electron occupation to 
toggle between and 1 electron. The energy of the trap 
is also fixed by the same gate voltage used to control the 
dot, but with a different capacitive coupling. This purely 
electrostatic model is a particular case of coupled quan- 
tum dots [9]. We chose arbitrarily to locate the trap on 
the source side. The source is grounded and a transport 
voltage is applied to the drain (Vd). We neglect any en- 
ergy dependence of the tunnel transparencies. The dots 
are neither simply in series, as there is a finite trans- 
parency between the source and the dot, nor in parallel 
as there is no direct source-drain current via the trap. 
The electrostatic model is fully characterized by the ca- 
pacitive couplings (see Fig. [1]): Cmn is the capacitance 
between m and n where m,n=d (dot), t (trap), s (source), 
dr (drain) or g (gate). Ct = Ctd + C^t + Cgt is the total 
capacitance of the trap. Cd = Ctd + Csd + C'drd + Cgd is 
the total capacitance of the dot. We consider sequential 
tunneling events driving the system into a finite number 
of states, defined by (Ndot, Ntrap), the occupation num- 
bers of the dot and trap. We deduce the drain-source 
current from the stationary occupation probability of the 
trap and the probability distribution of each charge state 
in the dot calculated with the rate equations. It is pos- 
sible to null the dot-trap or source-trap current with- 
out changing the conclusions. This does not change the 
electrostatic scheme but changes the resulting differential 
conductance lines (not discussed here). The reference of 
electrostatic energies is arbitrary and does not affect our 
model which implies only energy differences between con- 
figurations. 

Figure [2] shows the result for Cgt , Cgd , Ctd , Cst , Cd 



= 0.004, 16, 1, 2, 117 aF respectively i.e. for very small 
trap-gate capacitance. The current via the trap is 1 pA, 
to be compared with the dot-source and dot-drain con- 
ductances of 0.1 ^ (0.4 nA for Vd = lO""^ V). The level 

c 

arm parameter for the trap, at = — 0.0013 is very 
small, about 10 times smaller than in ref. [3]. The num- 
ber of successive distorted Coulomb diamonds increases 
as at decreases. A small at means that the trap is much 
less sensitive to the gate voltage than to the dot occu- 
pation and drain voltage. As we will see later the limit 
^ is similar to the DOS model. A small at also 
implies that the differential conductance lines are paral- 
lel to the diamond edges. We have checked that larger 
at gives different slopes. We set to volt the degeneracy 
point where the trap mean occupation number is ^ . Near 
this point we observe distorted Coulomb diamonds: two 
replicas appear, corresponding to the two charge states 
of the trap. The shift in gate voltage for the replicated 
diamond is given by AVg = q^^q^ • The largest distortion 
AV^ :^ occurs when Cgt = Ctd [3]. This occurs if 
the trap has equal couplings to source and dot. If the 
trap is more coupled either to the source or to the dot 
the apparent shift is reduced. 

A very important point is that the drain current van- 
ishes at small bias near the degeneracy point, when the 
diamonds are replicated. Figure [3] illustrates the physical 
origin of the effect: Starting from = V at small bias 
{Vd — V) when the gate voltage increases an electron is 
transferred from the trap onto the main dot. This cor- 
responds to the (N-1,1)^(N,0) transition (left square in 
figure [3]), which occurs before the (N-1,1)^(N,1) transi- 
tion (circle on figure [3]) because of the repulsion by the 
charged trap. This charge exchange between the trap 
and the dot does not give any drain-source current. More 
complex sequences with the same initial and final states 
and an intermediate state could in principle give a fi- 
nite current. First, co-tunneling events are not taken 
into account. One could also expect an extra electron 
to tunnel first from the source onto the dot, then the 
trap releasing its electron to the source. Alternatively 
the latter could happen first, then an electron could tun- 
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FIG. 3: Same simulation than in figure [21 focusing on tfie 
degeneracy region of the trap {Vg OV). Top panel: energy 
of lowest states at Vd = OV. Broken and solid lines are re- 
spectively for empty and occupied trap. Bottom panel: corre- 
sponding stability diagram. When increasing the gate voltage 
the (N-1,1)^(N,0) transition (left square) is reached before 
the (N-1,1)^(N,1) transition (circle). The system then fol- 
lows adiabatically the lowest energy state and the next transi- 
tion is (N,0)^(N,1) (right square). The (N-1,1)^(N,1) tran- 
sition is then replaced by (N-1,1)^(N,0)^(N,1). The latter 
does not give a drain current. 



nel from source to dot. However both sequences involve 
intermediate states too high in energy ((N,l) and (N- 

1.0) respectively), and therefore are forbidden. For the 
same reason, once in the (N,0) state it is impossible for 
an electron to exit the dot into the drain, because the 
(N-1,0) state is too high in energy at this gate voltage. 
Increasing the gate voltage further, the system reaches 
the (N,0)^(N,1) degeneracy (right square in Figure [3|). 
At this point the source can release an electron to the 
trap. This second transition at constant number of elec- 
trons in the dot does not yield any source-drain current. 
In summary the (N-1,1)^(N,1) transition has been re- 
placed by the (N-1,1)^(N,0)^(N,1) sequence where the 
trap occupation oscillates in gate voltage by successive 
transfers from the source and into the dot. In this gate 
voltage range the only way to recover a drain current is 
by applying a sufficient bias to allow N-1 or N electrons 
on the dot. Although a comparable situation has been 
described qualitatively in Ref. our calculations do 
not predict the 'kinks' they observe, but replicas instead. 

Far from the trap degeneracy (Vg OV), our model 
recovers undistorted Coulomb diamonds, as shown in 
Figure O Drain current is restored at low bias as (N- 

1.1) ^(N,l) is now the lowest energy transition. More 
interestingly, our model predicts lines of differential con- 
ductance in the non blockaded regions, both negative and 
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FIG. 4: Same simulation than in figure [21 focusing far from 
the degeneracy region of the trap [Vg ^ V). Here the trap is 
always occupied at small bias voltage and the diamond is not 
distorted. Top panel: energy of lowest states at Vd — 0.6 mV. 
Bottom panel: corresponding stability diagram. Region a : 
transitions (N-1,1)^(N,1) can occur. Both states are sequen- 
tially obtained, hence a current flows through the dot. Region 
c: in addition, (N-1,1)^(N,0) is also possible. When this oc- 
curs, no current passes through the dot anymore until the slow 
rate (N,0)^(N,1) transition happen. The state of the trap is 
a "dark state" lowering current through the dot. As a result a 
negative differential conductance line separates regions a and 
c. Region h\ like in region c the trap can exchange its electron 
with the dot, in addition the (N-1,0)^(N,0) transitions can 
happen. Therefore, whatever the trap state current through 
the dot can circulate and there is more current through the 
main dot than in a, hence a positive differential conductance 
line separates regions a and h. Region d: the dot can even be 
filled with N+1 electrons when the trap is empty. 



positive, as readily visible in Figure [2] and shown in more 
detail in Figure [H At some bias, the mean occupation 
of the trap is allowed to vary. The system can there- 
fore be in various charge states, implying different total 
conductances, hence maxima of differential conductance. 
We now discuss in more detail the origin of a negative 
differential conductance line between regions labeled a 
and c in Figure [H As mentionned before, in both of 
these regions (N-1,1)^(N,1) transitions are allowed and 
responsible for a finite drain current. In addition, the 
(N-1,1)^(N,0) transition is possible only in region c, al- 
lowing the trap to transfer its electron to the dot. When- 
ever this event happens no current passes through the dot 
anymore until the trap is filled again by an electron from 
the source. The latter event is much slower than tun- 
neling between source, dot and drain, so the current is 
smaller in region c than in a and therefore the differential 
conductance is negative. The (N,0) state (trap empty) is 
a "dark state" blocking the current through the dot. We 
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have checked that slowing the (N,0)^(N,1) transition by 
decreasing the source-trap tunnehng rate reduces further 
the current in region c. 

The main conclusion of this study over a large gate 
voltage range (Figure [2]) is that both the negative differ- 
ential conductance lines far from degeneracy and the sup- 
pression of current at low bias near degeneracy have the 
same origin, namely the electrostatic interaction between 
the trap and the dot. This effect can be exploited to build 
Coulomb blockade rectifiers suggested in ref. [HI, [lH . 



III. COMPARISON WITH THE ES AND DOS 
MODELS 

Usual explanations for lines parallel to the edges of 
Coulomb diamonds are based on the ES or DOS models. 
Transport excitation spectra have been investigated in 
much detail [13, [3] and their modification due to the 
Zeeman effect |3, photons [HI or phonons [16j are widely 
observed in the Coulomb blockade stability diagram of 
quantum dots. 

The ES Model is based on the existence of a resolved 
excitation spectrum because of the quantized kinetic en- 
ergy of electrons confined in a small volume. It requires a 
temperature much lower than the mean energy level spac- 
ing A. We illustrate how lines appear in this model in 
figure O The left panel is the orthodox model simulation 
of a quantum dot with A <C ksT. The charging energy 
is constant and equal to 1.38 meV {Cd = 116 aF). The 
right panel corresponds to a dot with the same charg- 
ing energy, at the same temperature but with 4 electrons 
occupying 4 non degenerate orbital states with spacings 
of 0.4, 0.3 and 0.1 meV respectively. Positive differen- 
tial conductance lines appear. In this constant charging 
energy model the one-particle excitation spectrum shifts 
between successive diamonds as levels get occupied, the 
first excited state for N electrons becoming the ground 
state for N+1 electrons. Very often scrambling of the 
one-particle spectrum after adding electrons is observed 
but recently a clear correlation between successive dia- 
monds has been reported [17]. The absence of scrambling 
is interpreted as the absence of variation of the dot shape 
when adding electrons. That requires a steep enough con- 
finement potential, for instance sharp etched edges. In 
the ES model it is natural to expect only positive differ- 
ential conductance lines because excited states (at higher 
energies than ground states by principle) are more tun- 
nel coupled to the electrodes. If equal tunneling rates are 
supposed for the ground and excited states, as done in 
Figure O then differential conductance is always positive. 
To observe negative differential conductance within the 
ES model the excited state must be less conducting. This 
could occur because of a specific fluctuation of the wave- 
function envelop, a selection rule [18], a blocking state 
[3 or Stark effect 

The DOS model is appropriate if the very same pat- 
tern is observed in successive diamonds [U, [22[ , in strong 
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FIG. 5: Left panel: Simulation of the stability diagram (differ- 
ential conductance versus gate and drain voltage) for a quan- 
tum dot, where the mean level spacing is negligible compared 
to the temperature T= O.IK (orthodox model). The charging 
energy is constant and equal to 1.38 meV (Cd= 116 aF). Right 
panel: simulation of a dot with the same charging energy at 
the same temperature, but with 4 electrons occupying 4 non 
degenerate orbital states with spacings of 0.4, 0.3 and 0.1 
meV respectively (constant charging energy model). Positive 
differential conductance lines appear outside the Coulomb di- 
amonds. The one-particle excitation spectrum shifts between 
successive diamonds as levels get occupied, the excited state 
for N electrons becoming the ground state for N+1 electrons. 



contrast with the ES model, because the energy in the 
electrodes does not depend on the gate voltage. In this 
model the source-drain current is proportionnal to the lo- 
cal density of states vs,d{E) in the source and drain. If 
for simplicity we suppose that the conductance is domi- 
nated by source-dot tunnel barrier, then G{V) ex In 
that case lines of differential conductance are expected to 
be parallel to the negative slope edge of the Coulomb dia- 
mond. This edge corresponds to the alignment of the last 
unoccupied level in the dot with the Fermi level in the 
source. In the opposite case of a dominating dot-drain 
barrier the lines will be parallel to the positive slope and 
G{V) oc ^. 

In heavily doped semiconducting electrodes, just like 
in any diffusive reservoir, the local DOS fluctuates with 
the energy because of quantum interferences of elasti- 
cally scattered quasi-particles diffusing coherently within 
a length scale related to their lifetime at a particular en- 
ergy |22|. The characteristic energy for these fluctua- 
tions is set by the inverse of the quasi-particle relaxation 
time, which decreases as the energy moves away from the 
Fermi energy. The local DOS fluctuations increase when 
the diffusion coefficient D = \v\t decreases, with d the 
dimensionality of the electrodes, vp the Fermi velocity 
and r the elastic mean free time. Indeed a very large D 
corresponds to a near perfect electronic reservoir without 
local fluctuations of the DOS. The local DOS fluctuations 
are reinforced in small, confined geometries. If the trap 
considered in our model gets closer to the source, at some 
point electronic orbitals in the source and trap will over- 
lap and a strong local fluctuation of the DOS appears due 
to hybridization. Then the energy of the trap depends 
less and less on the gate voltage. Therefore there is a con- 
tinuous evolution from the DOS model (which is a purely 
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quantum mechanical model) to our electrostatic model 
when increasing the trap-electrode distance. The DOS 
model can explain other features which cannot be under- 
stood within the ES model [23]. First the DOS model 
naturaly explains negative differential conductance lines, 
as the sign of changes. In average one expects as 
many positive and negative differential conductance lines. 
Also the DOS model, can explain lines pointing to ener- 
gies below the first available level in the dot. Such lines 
cannot arise from the ES model, as shown in the first 
diamond in the right panel of Figure [51 

Our model differs from the ES and DOS models by the 
3 following facts. 1) Lines are shifted progressively from 
one diamond to the next, in contrast to both ES and 
DOS models. 2) Lines coexist with sawtooth distortions 
(replicas) of diamonds at different gate voltages. 3) Neg- 
ative differential conductance appears at finite bias volt- 
age together with anomalously small current at low bias 
in distorted diamonds. Like the DOS model it explains 
why lines can appear at energies extrapolated below the 
ground state of the artificial atom. Finally our model, 
unlike both ES and DOS models, can explain lines not 
parallel to the diamond edges (if at is large) [J, |2j] . 



IV. SILICON QUANTUM DOT WITH ARSENIC 
DONOR AS A TRAP 

We originally developped our model for silicon 
nanowire transistors with implanted arsenic traps in the 
tunnel barriers [3]. Here we report new data recorded 
in similar samples which clearly show the correlations in 
successive undistorted Coulomb diamonds predicted by 
our model. Fig. [6)3 shows a typical stability diagram with 
a pattern of differential conductance lines very weakly de- 
pendent on the number of carriers in the dot. This weakly 
shifted pattern is a characteristic signature of our model 
for traps much more coupled to the source (or drain or 
dot) than to the gate. The samples are described in ref. 
[l|, [3, [Hi and schematically drawn in Fig. [6^. They are 
SOI-MOSFETs adapted in terms of doping to become 
controlled SETs at low temperature. A 20 to 80 nm wide 
wire is etched to form the channel. The source and drain 
parts of the wire are highly doped to form good metallic 
reservoirs (As, 0^ 10^^ cm~^). The central part of the wire 
is covered by a 40 nm poly- Si gate electrode, isolated by 
Si02, and self- aligned silicon nitride spacers (35 nm long) 
are deposited on both sides of the gate (rounded walls in 
Fig. [6^). The part of the wire below the gate and spacers 
(light gray regions on Fig. [6^1) is only lightly doped (As, 
5 X lO^^cm"^), so that at low temperature it forms an 
insulator. However directly below the gate electrode it 
can be tuned into a metallic state by applying a positive 
gate voltage. That way a quantum dot is formed under 
the gate. The tunnel barriers are the low-doped parts of 
the wire adjacent to the dot. The arsenic dopants inside 
these tunnel barriers are weakly capacitively coupled to 
the gate and are the traps considered in our model. Well 
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FIG. 6: (a) Schematic view of our MOS-SET. (b) Color plot 
of the drain differential conductance versus gate and drain 
voltages at T=350mK, which exhibits the same lines of con- 
ductance in successive diamonds with a small shift from di- 
amond to diamond. These lines are explained with a purely 
electrostatic model involving a background charge, (c) Simu- 
lation with our model and parameters given in the text. We 
obtain a line at the same position than in the experimental 
data and which evolves slowly with Vg because the trap is 
weakly coupled to the gate. 



centered donors give replica for the diamonds [3]. There 
is a gradient of Arsenic concentration (typically one order 
of magnitude for a 5 nm lateral distance in our process) at 
the border between the source-drain and channel. Many 
arsenic donors are located close, and are therefore well 
coupled to the source or drain. Such donor have a small 
lever arm parameter at and give undistorted diamonds 
with lines of differential conductance. Donors which are 
strongly coupled to the dot (due to electrostatic bending 
of the impurity band below the gate edges) produce the 
same effect. 

We performed our measurements at 350 mK in a ^He 
refrigerator and measured the differential conductance 
with a standard ac lock-in technique. At this temper- 
ature we do not expect to resolve the quantum lev- 
els in the dot. The mean energy level spacing A be- 
tween quantum states is the largest for small dots at 
low gate voltages where only a 2D electron gas is formed 
at the surface of the channel. In this limit we expect 
A2D ~ ~ 150 /ieV, with d = 4 the spin and val- 

ley degeneracy, m* = 0.19 mg the 2D effective mass, and 
A 4000 nm^ the total surface area of the gate/wire 
overlap, including the flanks of the wire. As the dot gets 



6 



filled, the electron gas eventually fills up the whole vol- 
ume of the wire below the gate and A falls below 20 jaeV. 
Quantum levels can only be resolved when A is larger 
than the width of the resulting lines of differential con- 
ductance. These lines have a full width at half maxi- 
mum of approximately 3.5 /cbT ~ 100 fieV given by the 
Fermi distribution in the leads. For the large gate volt- 
ages shown in figure [6] we are in the high density regime 
where A is too small to play a significant role. There- 
fore the sharp lines of differential conductance seen on 
figure [6] cannot be explained within the ES model. The 
lines have a typical energy separation of 1 meV, much 
larger than the calculated mean spacing of 20 /ieV. They 
are also observed at larger energy than expected in the 
DOS model. In our samples the reservoirs are highly 
doped silicon wires in which the mean level spacing is 
very small, and local DOS fluctuations have a correla- 
tion in energy of the order of the inverse inelastic time 
1 ns (— 4/ieV) at T=1K in heavily doped silicon 

In summary, unlike the ES and DOS models, our model 
explains quantitatively the weakly shifting lines between 
successive diamonds measured in Figure [6)3, with a trap 
located on the drain side of the dot. Figure [Bt shows 
the result of the simulation with Cgt, Cgd, Ctd, C'drt, 
Cd = 0.006, 13.3, 0.4, 0.046, 53.3 aF. 

V. CONCLUSIONS 

We have extended a simple electrostatic model of a 
charge trap coupled to a dot to the case of very weak 
coupling to the gate. In this regime new features are 
predicted over a large gate voltage range. Near degen- 
eracy of the trap the sawtooth pattern calculated in a 
previous work is recovered and the current suppression 
at low bias voltage is understood in more detail. We 



obtained new features far from this degeneracy, where 
Coulomb diamonds are not distorted. Lines of differen- 
tial conductance appear in the diamonds, very similar 
to the ones usually attributed excited states, although 
our model does not involve a discrete spectrum for the 
dot. These differential conductance lines can be posi- 
tive or negative, and parallel to either edge of the dia- 
mond. They also evolve very weakly with gate voltage, 
an original feature not predicted by other models. Our 
model easily accounts for negative differential conduc- 
tance lines, a feature usually attributed to density-of- 
state fluctuations in the contacts. This model and ours 
are converging when considering a trap located very close 
to the electrodes, as hybridization occurs and coupling 
to the gate goes to zero. Although the most basic signa- 
ture of our model is the sawtooth pattern and associated 
current suppression, we emphasize that its experimental 
observation is not required to validate our charge trap 
scheme. Indeed lines can very well be observed in undis- 
torted diamonds while the degeneracy region remains out 
of the energy range which can be probed experimentally. 

Even though our model involves only a single trap oc- 
cupied with zero or one electron and a dot treated in the 
orthodox model, it already gives a complicated pattern 
of lines and features in the stability diagram. It provides 
a quantitative but simple way to simulate the Coulomb 
blockade spectroscopy of quantum dots, and shows the 
great impact of a single charge on this spectrum. Further 
extensions like several traps, resolved mean one-particle 
level spacing in the dot, no n- negligible current through 
the trap, double occupation of the trap, Zeeman effect on 
the trap energy can be implemented in the near future. 
As more and more nanostructures designed for transport 
experiment exhibit Coulomb blockade, our model could 
account for many features observed experimentally as the 
presence of charge traps is very realistic. 
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